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Abstract. The calculation of the correlation functions of Bethe ansatz solvable models is 
very difficult problem. Among these solvable models spin 1/2 XXX chain has been investi- 
gated for a long time. Even for this model only the nearest neighbor and the second neighbor 
correlations were known. In 1990's multiple integral formula for the general correlations is 
derived. But the integration of this formula is also very difficult problem. Recently these 
integrals are decomposed to products of one dimensional integrals and at zero tempera- 
ture, zero magnetic field and isotropic case, correlation functions are expressed by log 2 and 
Riemann's zeta functions with odd integer argument C{3), C(5), C(7), .... We can calculate 
density sub-matrix of successive seven sites. Entanglement entropy of seven sites is calcu- 
lated. These methods can be extended to XXZ chain up to n = 4. Correlation functions are 
expressed by the generalized zeta functions. Several years ago I derived new thermodynamic 
Bethe ansatz equation for XXZ chain. This is quite different with Yang- Yang type TBA 
equations and contains only one unknown function. This equation is very useful to get the 
high temperature expansion. In this paper we get the analytic solution of this equation at 
A = 0. 
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1 Introduction 

We consider the spin 1/2 XXZ chain 



1=1 ^ ^ 1=1 



(1) 



Among the solvable models this model has been investigated for a long time. 

It was believed that the exact calculation of correlation function is impossible except the 
nearest neighbor correlation function. I derived the second neighbor correlation function for 
J<0, T = 0, /i = 0, A = l using the Lieb-Wu solution of one-dimensional Hubbard model |1]. 
Details are given in Appendix [A) In 1990's multiple integral formula were proposed for h = 0, 
T = and recently multiple integral formula was extended to h ^ 0, T ^ 0. But the factorization 
of these multiple integrals to the integrals of lower dimension still remains a difficult problem. 
I explain the present situation of factorization in Section [2l 

Since Yang and Yang proposed the thermodynamic Bethe ansatz (TBA) equation for one 
dimensional bosons [2], the calculation of free energy becomes important for other solvable 
models. Yang- Yang type integral equations for ([1]) were proposed in early 1970's [3 HI [5]. In 
this theory infinite number of unknown functions appeared for |A| > 1 and finite number of 
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unknowns appear for A = cos ttu with u = rational number. The equations change by the value 
of A and not so convenient for numerical calculations. But some important physical properties 
at low temperature were investigated using these equations. 

Around 1990 the quantum transfer matrix method was applied to this model and numerical 
results coincide with those of Yang- Yang type equations [3 O [8]. In 2001 I proposed a new 
TBA equation which contains only one unknown function for ([T]) [9l |TOj . This equation is very 
convenient to do the high-temperature expansions and we get one hundred-th order of high 
temperature expansion [11]. The traditional cluster expansion method gives up to 22'nd. This 
method was also applied to the Perk-Schultz model [12]. 

Unfortunately this equation does not numerically converge at low temperature like T/\J\ < 
0.07 because the integrand strongly oscillates. Some other numerical method is necessary. 

1 showed that this equation gives the known exact results in Ising limit (A — )■ cxd and J A = 
finite) [9]. In Section [3] I give the analytic solution of this equation for XY case (A = 0). 

2 Correlation functions of XXZ chain 

2.1 Known exact results for J < 0, T = 0, h = 0, A = 1 

1. Nearest-neighbor correlator 

(5|5|, i) = i- - iln2 = -0.1477157268. . . (2) 
12 3 

from the ground state energy per site by Hulthen |13j (1938). 

2. Next nearest-neighbor correlator for XXX 

(5|S|+2) = - ^ in 2 + ^C(3) = 0.06067976995 . . . , (3) 

from the ground state energy of the half-filled Hubbard model by Takahashi [1] (1977). 

3. The twisted four-body correlation function 

{{Sj X Sj+i) ■ {Sj+2 X Sj+3)) = ^ In 2 - ^C(3) = -0.104197748 (4) 
from the third derivative of transfer matrix by Muramoto and Takahashi |14j (1999). 

2.2 Multiple-integral representations for density matrix 

For more general XXZ model with an anisotropy parameter A 

1. A > 1, T = 0, /i = 0: Vertex operator approach Uq{sl{2)) Jimbo, Miki, Miwa, Nakayashiki 
[111 US] (1992). 

2. -1 < A < 1, T = 0, /i = 0: qKZ equation Jimbo and Miwa [IT] (1996). 

3. Rederivation by the quantum inverse scattering method. Generalization to the XXZ model 
with a magnetic field, T = Kitanine, Maillet, Terras [IB] (1998). 

4. Integral formula for finite temperature and finite magnetic field Gohmann, Kliimper, 
Seel [E] (2004). 

Correlation function for successive elementary block 
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0/,,' J \o 0/,,' J \i o;,,' J Vo i/ri' 

is represented by n-fold integral. For example the emptiness formation probability (EFP) for 
XXX chain at T = /i = is 



(-vr)^2-"/ d-Ar[ "7^'^";~"''m 



A- sinh7r(Aa - Ab) -A- (A^ - ^ {Xj + ^ 



Other arbitrary correlation functions over successive n-sites have similar n-fold integral repre- 
sentation. 

2.3 Boos— Korepin method to evaluate integrals 

Here I introduce the details of direct factorization of multiple integrals by Boos and Kore- 
pin [201 El]. 

1. Transform the integrand to a certain canonical form without changing the integral value. 

2. Perform the integration using the residue theorem. 

Example: 

3 .oo-i 



/°°~ 2 HA ■ 
. y^f/3(Al,A2,A3)r3(Ai,A2,A3), 

n sinh7r(Aj - Afc) 



TsCAi, A2, A3) — 



3 ' 
Y\ sinh^vrAj 

(Ai + i)2A2(A2 + i)Ai 



(A2 - Ai - i)(A3 - Ai - i)(A3 - A2 - i) ' 



Details of transformation are given in Appendix [Bl 
In canonical form denominator is 

n(^2fc-i -A2fc), l<[n/2]. 

k=l 

The n-dimensional integral is decomposed to one and two dimensional integrals. Transform 
T3(Ai, A2, A3) into a canonical form Tg (Ai, A2, A3) 

^ , p(3) 

T3(Ai, A2, A3) ~ r3=(Ai, A2, A3) = + ' 



A2 — Ai ' 



pP = -2X2X1 Pf ^ = --iXi- iXs - 2A1A3. 
Perform the integration 

4' 



1 
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P(3) = Jf + Jf = i-ln2 + ^C(3). 

Thus second neighbor correlator was rederived from the integral formula. 
For P(4) = p^\^\ the integrand is 

(A43 - 0(-^42 - 0(-^41 - i)(A32 - «)(A31 - «)(A21 - i) ' 

Here we put A^fe = A^ — A;,. Canonical form is Pq^^ + P^^V^^si + P2^V('^2iA43), 

p(4) _ 34 2x3 

p[^^ = A?(30A|A^ + 30iA3A| - 16A^ + ISAsA^ + 8A4) 

+ Ai(30iA|A^ + 30A3A^ - 16iA^ + ISiAgA^ - A\l + 4^4) 

- 20AiAl - 20a3Al + ^Al - ^AgAf - ^iA4, 

5 5 lU 

Pi") = 2A?Ai + 4UiAi - ^Ai - ^AiAs - iAg + \, 

P(4) = 1 - 21n 2 + — C(3) - - In 2 • C(3) - -C'(3) - -C(5) + -ln2 • C(5). 

In 2003, we calculated p^j^I^I by Boos-Korepin method and obtained all the correlation func- 
tions on 4 lattice sites [22]. Especially, the third-neighbor correlator is 

<5|S|+3> = 1 - 31n2 + ^C(3) " ^1^2 • C(3) - ^C(3)^ - ^C(5) + f 1^2 • C(5) 
= -0.05024862725... . 

The other correlation functions for n = 4 are 

11 qi 1 

<5JSjVi5|+25|+3> = ^ + 1^1-2 - — C(3) + -ln2 . C(3) 

+ 1^(3)2 + ^^(5) _lln2.C(5), 
80^^ ' 96^^ ' 24 

(S^S^,.S^,nSL.) = ^ - lln2 + — C(3) - — ln2 • C(3) 
\ J j+i j+2 j+3/ 240 6 120 ^12 ^ ' 

_ 1^(3)2 _^<^(5) + 5^n2. C(5), 
20^ ^ 96^ ^6 ^ ^ 

11 1 69 5 

(5|S|+,5|+25jV3> = ^ - + — C(3) - -ln2 . C(3) 

_ 1^(3)2 _^<^(5) + 5^n2. C(5), 
20^ ^ 96^ ^6 ^ ^ 

/ QZ QZ QZ \ / QX QX QZ QZ \ _|_ / QX QX QZ \ I / QX QZ QZ QX \ 

\ i i+1 j+2'-'j+3/ — \'^j'^j+l'^j+2'^j+3/ "T Wj '^j+l'-'j+2'^j+3/ "T X'-'j '-'i+l'-'j+2'-'j+3/ ■ 

From these results we can reproduce the twisted correlation function in Q. 
In a similar way, P(5) was calculated after very tedious calculations [23] 

But the direct integral of other correlations for five sites is almost impossible. 
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2.4 Algebraic approach and qKZ relation 

Next problem is to calculate {SjSj_^^) for XXX model. In principle, it's possible to calculate 
other five-dimensional integrals by use of Boos-Korepin method. It, however, will take tremen- 
dous amount of time. 

We propose a different method ("algebraic approach") and obtain analytical form of (SjSj^^). 
This is a generalization of the method by Boos, Korepin, Smirnov (2003) for P(6) [23]. We 
consider the density matrix p of successive n sites of inhomogeneous six vertex model with 
different spectral parameter Zj for j-th site 

lim Pn\ei%Jzi,Z2, ■ ■ ■ , Zn) = {EII ■ ■ ■ Et^) , 
{Ei)s^s> = Se,sSe',s', €, e' = ±1. 



For n = l,and 2 we have 

~ 2' 



Pt, + -{Z1,Z2) = -^UJ{ZI -Z2), P2^ + _(2l,Z2) = ^^{Zl - Z2), 

with 

u^i.) ^ \ + 2Y,{-lfk^^ = i - 2(1 - x^) ^ x^\.{2k + 1) 

k=l fc=0 

00 

Ux) ^ J](-l)"-V-- = (1 - 2i--)C(x), Ca(l) = ln2, 

n=l 

x{x + 2) , ^ 3 1 

— 7, Ujix) 

x2 - 1 ^ ' 21 -X 



uj{x + l) = ^ —uoix)--- ^, uj{-x) = uoix), a;(±zoo) = 0. 



The general element of density matrix must satisfy the following algebraic relations. 

• Translational invariance 

Pn,ti,...,fin{zi -\- X, . . . , Zn -\- x) = Pn,ei,...,€n{zi, ■ ■ ■ , Zn}- 

• Transposition, negating and reverse-order relations 

Pn,ei,...,€„{zi, ■ ■ ■ , Zn) = P„'i ' ,/ (~2;i, . . . , — Z„,) 

— e' — e' e' e' 

~ Pn,-ei,...,-e„(^li • • • ' ^n) = Pn',t„,..},ei{ — Zn, • • • , — ^l)- 

• Intertwining relation 

R-j-^,^^{Zj - Zj + l)p,„e^ll^e]...{- ■ ■ Zj + l,Zj . . .) 

j j + l 



i?++(^) = i?Zl(^) = 1, RlZ{z) = RzXiz) = 
Rt^iz) = R~tiz) = 
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• Reduction relation 

Pn,+'y2','.'."E„(^l> ^2, ■■■,Zn)+ Pnf- 'e2,%ni^l,Z2, • • • , Zn) = Pn-Tll...,eA^2, ■ ■ ■ , ^n)- 

• First recurrent relation 



5 • • • 5 -^nj} 



Pn\fue2,%Az -l,Z,Z3,...,Zn) = -^e\,-t'^^\^2Pn-tlZ^^^..,eS^^ ^3 

Second recurrent relation 

lim Pn';l'y2','.'"en(^l, ^2, ■■■■,Zn) = b J ]-p^^::'{l^:_ {Z2, . . . , Zn). 

Identity relations 



, . . . , Zyi). 



^ ^ Pn,ei,.."e„(^l) • 



^ ^ /5n,ei,.."en(^l) • • • ) ^r! 



£1, •••,<;« 



= Pn';+,t.,+ {zi, ..■,Zn)= Pn;-1..-{ZI, ■ ■ ■ , ^^n)- 

If we assume as follows 

P3^+++{zi,Z2,Z3) = +^(zi,Z2k3)'^(2l - Z2) + A{zi,Z3\z2)uj{zi - Z3) 

+ A{Z2,Z3\Z2)UJ{Z2 - Z3), 
P3~-++ (^1 , ^2 , 23 ) = P2^+ (^2 , ^^3 ) - /53;^++ (^1,^2,23), 
PZ~-+-{Z1^Z2,Z3) = P2~^+{zi, Z2) - P3~^++{zi, Z2, Z3), 
{Zi - Z3)iz2 - Z3) - 1 



A{zi, Z2\Z3) 



12(zi - Z3){Z2 - Z3) 



these relations are satisfied. In the homogeneous limit Zj — >• this gives the correct correlation 
functions of XXX model. In the homogeneous limit each term diverges but we have finite 
limiting number. Then we can calculate the correlation functions of arbitrary element of density 
matrix using these algebraic relations, although the calculation become complicated. We have 
calculated all the inhomogeneous correlation functions up to n < 4 from the multiple integrals 
and confirmed these relations are fulfilled. 

Further we have found the inhomogeneous correlation functions can be represented in terms 
of w-function 



\j = l / "1=1 I<fcl<fe3<fc5<---<fc2m-l<n, k2m>k2m^l 

e' e' 

AeY;::;,el{ki, . . .,k2m\zi, . . .,Zn)uj{zk-, - Zk^) ■ ■ ■ Uj{zk2rr,-i - Zk2m)^ 
^ei,...,e„ (^1 1 • • • ) k2m\zi , . . . , Z^) 



£■,,...,€„ |_ _ Qei,...,en{kli ■ ■ ■ , k2m\zi, ■ ■ ■ 1 Zn) 



: rational function of zi, . . . , Zn- 
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Denominator is 

2m / \ 

n (^^2,-1 - ^fc2i-l)(^fc2,-l - ^fc2i)(^fc2, - ^fc2;-l)(^fc2, " ^^2; ) H I ll ~ I " 

l<j<k<m 1=1 \ ij^ki,k2,...,k2m / 

The total exponent for this is 4nm — 2m? — 2m. The largest exponent for Zi is n — 2 for 
i G {ki, . . . , k2m} and 2m for i ^ kj. Numerator is also polynomials of zi, . . . , Zn which satisfies 
the same exponent conditions. Unknowns are the coefficients of each terms. Number increases 
drastically as n, m increases. Algebraic relations give the over complete linear equations. By 
using mathematica we have unique solution of these equations. By use of algebraic relations, 
we have calculated all the polynomials for n = 5 

Qlu.'.'.'Mikl,k2\zi, . . . ,Z5), Qlu.'.T4ih,- ■ -^k^lzi,.. .,Z5). 

By the memory problem this calculation stopped at n = 6. 
Fourth-neighbor correlation function |25] 

- f C(3) . C(5) - HpC(5)H 5H5C(T) - l^ln 2 . C(7) + 5f C(3) . C(7) 
= 0.034652776982 . . . , 

(5|5|+5) = 3^ - f Ca(l) + ^Ca(3) - ^Ca(l)Ca(3) - ^(^(3)^ - 608^(3)^ 

- ^Ca(5) + ^Ca(l)Ca(5) - ^Ca(3)Ca(5) + 3648Ca(l)Ca(3)Ca(5) 

- ^C.(3)^C.(5) - ^a(5)^ + ^C.(1)C.(5)^ + ^a(3)C.(5)^ 
+ ^C.(5)3 + ^C.(7) + 56952C.(3)C.(7) - ^C.(1)C.(7) 

- ^^Ca(l)Ca(3)Ca(7) - ^^U^fCa{7) + 7840Ca (5)Ca (7) 

- 8960Ca(3)Ca(5)Ca(7) - ^^Ca(7)2 ^ 3136OC, (l)Ca (7)^ - 686Ca(9) 

+ 18368Ca(l)Ca(9) - 53312Ca(3)Ca(9) + 35392Ca(l)Ca(3)Ca(9) 

+ 16128Ca(3)2Ca(9) + 38080Ca(5)Ca(9) - 53760Ca(l)Ca(5)Ca(9) 
= -0.03089036664760932 .... 

Using this algebraic method we can calculate all the element of density sub-matrix of successive 
6-sites. Longer system is quite difficult because of the memory and computing time problem. We 
can calculate 6-th neighbor and 7-th neighbor correlations using the generation function method. 
They are represented by long polynomials of CaS. Here we write only numerical results |26j 

{SjS]+Q) = 0.02444673832795890 . . . , {S-S^+j) = -0.0224982227633722 .... 
2.5 Calculation by continuous dimensions 

In a series of papers Boos, Jimbo, Miwa, Smirnov and Takeyama formulated these algebraic 
calculation by the trace of continuous dimension of auxiliary space [271 EEl EU [30] 

(pn)!;;:::;!: = (vac|(i^|j)i • • • (i^i:)n|vac). 
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Kiei, en,en, ■ ■ ■ , ei) = (-1)" (Pn)ri'X.r'"> 



n ^ 

Sn = Y[:;{\+)j\-)j- |-)il+)j)> hn = exp(0„)s„. 



2 

i=i 

Density sub-matrix in 2" dimensional space is mapped to a vector in 2^*^ dimensional space. 
Qn is an operator in this space. Monodoromy matrix is defined as follows: 

r„(A) = Lf (A - zi - 1) • • • 4°^ (A - z„ - 1)lW(A - z„) • • • lS°)(A - zi). 

Here /, H, E, F are d x d matrices 

hj = = (d+l - 2i)6ij, Eij = {i- l)5ij+i, Fij = {d - (5) 

These satisfy commutation relations [H,E] = —2E, [H,F] = 2F, [E,F] = —H. aj, , aj are 
Pauli operators in 2n spin space j = 1, . . . , ?i, 1, . . . , n. Especially at c? = 2 we have 

-G?)' "-^o')' -K°o)' -(oJ 

and 

/^(A) B{\)\ 

^r^y"-) - \C{\) D{X))- 

The trace of any monomial of /, E, H and is a polynomial of dimension d. One can calculate 
from the definition ([5]). For example, 

Tr^I = d, TidHH = {d^ - d)/3, TraH = 0. 



The dimension d is replaced hy fi — u. Following [30], the operator 0,^ is given by 
1 f dfi dv (/i — v)uj{^ — v) 



" 2 J 2TTi2-ni 



X Tr^_,r„ (^(/^)i^(^) + D{^,)A{v) - B{^,)C{v) - C{ti)B{v)), 



where the integration path should surround all Zj counter-clockwise. In the homogeneous limit 
1 f dfi dv {n — iy)uj{fi — v) 



Zj ^ becomes 



2 7 27ri 27ri (1 - (/i - vYY^l'^il - /i2)"zyn(l - i/2)n 

X Tr^_,r„ (^(/^)^(^) + D{^l)A{u) - i?(/i)C(i.) - C(/i)i?(z.)), 

and the calculation becomes very simple. By this formulation we could calculate all the elements 
of density sub-matrix at n = 7. Calculating the eigenvalues of matrix, we can calculate the von 
Neumann entropy (entanglement entropy) up to seven sites, 

S{n) = -ir Pn log2 Pn = - ^^^a log2 ^a-, 

a=l 

where oja are eigenvalues of density sub- matrix In Tabled 5 (n) is given up to n = 7. 
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Table 1. von Neumann entropy S{n) of a finite sub-chain of length 



5(1) 


5(2) 


5(3) 


5(4) 


1 


1.3758573262887466 


1.5824933209573855 


1.7247050949099274 



5(5) 5(6) 5(7) 

1.833704916848315 1.922358833819333 1.997129812895912 



2.6 Generalization to XXZ model 

In 2003, Kato, Shiroishi, Takahashi, Sakai have generalized the Boos-Korepin method to the 
XXZ models with an anisotropy parameter |A| < 1 for successive three sites [31]. For example 
P{n) is represented as follows: 

P(n) = {-.r^ r r ^ rr sinh(x. - X,) 

y_oo 27r y„oo 2tt smh [{xa - Xb - ivr) u) 

^ sinh""J' {{xj + f ) z/) siny-i {{xj - if) u 

X 



n 



cosh" X 

Here A = cos(7ri/). Similar integral representations for any arbitrary correlation function for 
successive three sites were calculated. 
Nearest-neighbor correlation functions 

Next nearest-neighbor correlation functions 
Here 

Cj := cosTTjV, Sj := sinvrjV, 

^ , , f°^~^ 1 coshzvx ^, , , f°°~^ 1 (5 coshi/x , , 

(^(j := / dx — : , C'Jj) := / dx — — ■. . 6 

J_oo-^ smhxsinh-^ra J-oo-^ smhx di/ sinh-^ z^x 

Replacing i'n/'^j can also get the correlation functions in the massive region A = coshry > 
1 [32]. Third neighbor correlations is also expressed by functions (^^ and C^, although the 
expression becomes more complicated [33] • In Figs. [1] and [21 the nearest neighbor, the second 
neighbor and the third neighbor correlations are shown as functions of A. 

3 Simplified thermodynamic Bethe ansatz equation 

Simplified thermodynamic Bethe ansatz equation at temperature T [9l [10] is 
u{x) — 2cosh(/i/r) - 



2 



„ , 2ttJ sin 6 
coth-(x -y- 2i)exp ( — — ai(y + i) 



„ , 27rJsin0 
+ coth -[x-y + 2i) exp ( — — oi(y - i) 



1 dv 

0, (7) 



u{y) 27ri 



10 



M. Takahashi 



0.2 










-\ /qxox\ 
^^^\'-'0 '-'2 1 




1 


2 3 4 








-0.2 







A 



Figure 1. The nearest-neighbor and the next nearest neighbor correlation functions for the XXZ chain. 

We calculated (S'|5'|+i), (5'|S'|+2) and (S'JSJ+j)- 



where 



' sm I 



ai[x) 



27r(cosh 6x — cos 6) 
Free energy per site is 
/ = -rinn(O). 

In this section we look for analytic solution for XY case A = 0, = 7r/2. For this case 
equation (j7]) becomes 



u{x) — 2 cosh(/i/T) ~ 'j> ~^ tanh -^{x — y)2 cos 



J \ 1 dy 
Tsinh7r?//2y u(y) 27ri 

Putting X = tanh7rx/4, Y = tanh7ry/4 and using u{—Y) = u{Y) we have 

J i-y2\ Y 



0. 



u{X) -2cosh{h/T) + 2(1 - X^) f cos 



dY 



2T Y y (1-X2y2)(i_y2)^(y) 27ri 
Consider the Fourier transform of following function 



0. (8) 



g{e) = ^ In (^4 Tcosh^ (^sme]+ smh\h/T) 







27r 

Ztt 



Assume that 



lnu{X) = ao + 2^ajX2j. 



(9) 
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-0.25 L 

Figure 2. The third-neighbor correlation functions for the XXZ chain. 



We can show that this satisfies ([8]). This series is convergent at \X\ < 1. 

oo 

ln'u(e^^) + ln'u(e"*^) = 2ao + 2 ^ aj cos{2j9) 

= 2g{e) = In ^4 ^cosh^ (^^sm9^ + smh."^ (h/T) 
Then we have 

u{X)u{l/X) = 2cos (^{X - 1/X)j +2 + 4sinh2(/i/r). (10) 



The function u{X) has zeros at X = ziz(3j^±, where f3j^± = aj^± + + and Cij^± 

J 

oo ^ ^ 

X2 VI 



'^^'^-{j - ^) ± One can show that |/3j,±| > 1. u{X) should not have zeros at |X| < 1. Then 



we have u{X) = const • H ~ ~ l^^)- condition u(l) = 2cosh(/i/T) we have 

i=i 

infinite product expansion of u{X) 



oo (1-41) (1-41) 

n(X) =2cosh(VT)n773%T773%T- (H) 
From (jlOp we have 

2 cos {J^{y - l/^")) /u{Y) = u{l/Y) - (2 + 4sinh2(/i/T))/n(y). 

oo 

Then we can write u{X) = 2 cosh(/i/T) + v{X'^){l — X"^), v{X) = ^ djX'j . Convergence radius 



of is infinite 

l.h.s. of dlD = u{X) - 2cosh{h/T) + {l - X^) j u{l/Y) — 



Y dY 



{1-X-) 



dY' 
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Thus we have proved that ([TT]) . ([9]) satisfies the equation ([8]). The free energy 

-Tlnti(O) = -Too 
coincides with the known result |34] . 

4 Summary 

For J<0, A = l,T = /i = Owe obtained the factorized form of density sub- matrix up to 
n = 7. The entanglement entropy for seven sites is new result of this paper. Up to six sites we 
published in [26] . 

The six-th neighbor and the seven-th neighbor correlations are calculated by the generating 
function method for J < 0, A = 1, T = /i = [26]. 

For arbitrary A, T = /i = we obtained the factorized form up to n = 4. Correlations are 
given by two transcendental functions CuU) C(j) with j = 1,3, 5, . . . defined by dH). For 
correlations of n > 5 the calculation becomes very tedious and no one has succeeded. 

For simplified TBA equation, we obtained the analytic solution for XY limit A = 0. Analytic 
solution for Ising limit was given in |9]. 



A Strong coupling expansion of the Hubbard model 

The Hubbard Hamiltonian is written as follows: 



Na 



<ij> <y 



If we treat the interaction term as main Hamiltonian and hopping therm as perturbation in the 
half-filled case 

the effective Hamiltonian becomes as follows: 



i<j Li<j i<k 

+ tijtjktkitiii5{aj ■ (Jk){ai ■ ai) + 5(cTi • aj){ak ■ cri) - 5(crj • o-fe)(crj • ai) 

i<j<l,i<k,kj^j,l 



ai ■ Uj - Uj ■ ak - Ok ■ (Ti - (Ji ■ Ui - Ui ■ ak - aj ■ ai + l) 



CTi ■ Uj = 4Si • Sj. 

For one-dimensional half-filled case the four spin term disappears and the effective Hamiltonian 
becomes 

^ YiAS, ■ S,+i - 1) + ^ ^{4(1 - 4S, • S,+i) + (4S, • S,+2 - 1)} + O . (12) 

1 f' ^ 
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On the other hand exact ground state energy per site is expanded as [351 136] 

Jo{u))Ji{u))doj 



e = -4t 



oo 







= -4\t 



uo[l + exp{2U'uj)] 

2 /I o\ 2 



(13) 



Comparing the first term of (I12p and (jl3p . we get nearest neighbor correlation ([2]). Prom the 
second term we get the second neighbor correlation 

B Transformation to canonical form in case of -P(3) 

(Ai+z)2(A2 + i)A2Ai 
(A3 - Ai - i){X3 - A2 - i)(A2 - Ai - i) ■ 

is decomposed to the following three terms, 

i(Ai + i)2(A2 + i)A2A| ^ i{Xi + i)^{X2 + i)X2X3 i{Xi + i)^{X2 + i)X2Xl 



(A3 - Ai - i)(A2 - Ai - i) (A3 - Ai - i)(A3 - A2 - (A3 - A2 - i)(A2 - Ai - z) ' 
Using the antisymmetry of C/3(Ai, A2, A3) the first and the second terms are simplified as follows: 

r \2x (Ai + i)^A3 

the iirst term ~ Ai A2 — ^ ^ , 

A2 - Ai - z 

the second term ~ A?A2 — . 

X2 — Xi — i 

The third term is transformed as follows 

,2, . {Xi + i)HX3 + i)' , .(Ai+i)^Ai ziX^+zfXl 
the third term ~ —X1X2 — i ^ ^ : h i- 



X2 - Xi -i X2 - Xi -i (A3 - A2 - ■i)(A2 - Xi -i) 
Then T3 is transformed as follows: 

T -A a2 _ ijXi+ifXl 

^ ' (A3-A2-i)(A2-Ai-i)- 

We should note that the pole at A3 = and Ai = —i of U3 is canceled by numerator of the 
second term. So we can change the integration path Ai — t- Ai — i and A3 — )■ A3 + i 

I3 ^ — A2A3 



(A3-A2)(A2-Ai)- 
Using the antisymmetry of U we have 



X x2 3AfA| + 3iAiA| + 3iA?A3-Ai-A3Ai-A? 
~ -A2A3 . 

A2 — Al 

Thus the denominator is drastically simplified. Using the following relations 

-fih)-(=!^^-liXl+^)]f{X3), (14) 



A2 — Al \ A2 — Al 3 
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we can reduce the power of numerator, 

o i — iXi — iX^ — 2A1A3 

A2 — M 

Thus we have obtained the canonical form for -P(3). The derivation of (jl4p is as follows: 
-/(As) ~ -T ^ :J[X3) = — ; r -. T T T /(As) 



A2 - Ai A2 - Ai - i \ A2 - Ai - i A2 - Ai - i 

A3 _ (A, + ,)3 (^^^x /A|^A^+^ (A^+^^ 

+ ^ ^ J(As) ~ — ;^ r — + ^ ^ /(A3J 



A2 — Ai — i A2 - Ai / V A2 — Ai — i A2 - A 
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